EHMT2 Inactivation in Pancreatic Epithelial Cells Shapes the Transcriptional Landscape and Inflammation Response of the Whole Pancreas.

The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that EHMT2 inactivation alters growth and immune gene expression networks, antagonizing KRAS-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of EHMT2 in maintaining a transcriptional landscape that protects organs from inflammation. Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from EHMT2 conditional knockout animals ( EHMT2 fl/fl ) targeted to the exocrine pancreatic epithelial cells ( Pdx1-Cre and P48 Cre/+ ), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the EHMT2 fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels. The EHMT2 fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional EHMT2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific EHMT2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the EHMT2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

Structural and Dynamic Analyses of Pathogenic Variants in PIK3R1 Reveal a Shared Mechanism Associated among Cancer, Undergrowth, and Overgrowth Syndromes.

The PI3K enzymes modify phospholipids to regulate cell growth and differentiation. Somatic variants in PI3K are recurrent in cancer and drive a proliferative phenotype. Somatic mosaicism of PIK3R1 and PIK3CA are associated with vascular anomalies and overgrowth syndromes. Germline PIK3R1 variants are associated with varying phenotypes, including immunodeficiency or facial dysmorphism with growth delay, lipoatrophy, and insulin resistance associated with SHORT syndrome. There has been limited study of the molecular mechanism to unify our understanding of how variants in PIK3R1 drive both undergrowth and overgrowth phenotypes. Thus, we compiled genomic variants from cancer and rare vascular anomalies and sought to interpret their effects using an unbiased physics-based simulation approach for the protein complex. We applied molecular dynamics simulations to mechanistically understand how genetic variants affect PIK3R1 and its interactions with PIK3CA. Notably, iSH2 genetic variants associated with undergrowth destabilize molecular interactions with the PIK3CA receptor binding domain in simulations, which is expected to decrease activity. On the other hand, overgrowth and cancer variants lead to loss of inhibitory interactions in simulations, which is expected to increase activity. We find that all disease variants display dysfunctions on either structural characteristics or intermolecular interaction energy. Thus, this comprehensive characterization of novel mosaic somatic variants associated with two opposing phenotypes has mechanistic importance and biomedical relevance and may aid in future therapeutic developments.

Cancer-associated polybromo-1 bromodomain 4 missense variants variably impact bromodomain ligand binding and cell growth suppression.

The polybromo, brahma-related gene 1-associated factors (PBAF) chromatin remodeling complex subunit polybromo-1 (PBRM1) contains six bromodomains that recognize and bind acetylated lysine residues on histone tails and other nuclear proteins. PBRM1 bromodomains thus provide a link between epigenetic posttranslational modifications and PBAF modulation of chromatin accessibility and transcription. As a putative tumor suppressor in several cancers, PBRM1 protein expression is often abrogated by truncations and deletions. However, ∼33% of PBRM1 mutations in cancer are missense and cluster within its bromodomains. Such mutations may generate full-length PBRM1 variant proteins with undetermined structural and functional characteristics. Here, we employed computational, biophysical, and cellular assays to interrogate the effects of PBRM1 bromodomain missense variants on bromodomain stability and function. Since mutations in the fourth bromodomain of PBRM1 (PBRM1-BD4) comprise nearly 20% of all cancer-associated PBRM1 missense mutations, we focused our analysis on PBRM1-BD4 missense protein variants. Selecting 16 potentially deleterious PBRM1-BD4 missense protein variants for further study based on high residue mutational frequency and/or conservation, we show that cancer-associated PBRM1-BD4 missense variants exhibit varied bromodomain stability and ability to bind acetylated histones. Our results demonstrate the effectiveness of identifying the unique impacts of individual PBRM1-BD4 missense variants on protein structure and function, based on affected residue location within the bromodomain. This knowledge provides a foundation for drawing correlations between specific cancer-associated PBRM1 missense variants and distinct alterations in PBRM1 function, informing future cancer personalized medicine approaches.

Enhancing prognostic power in multiple myeloma using a plasma cell signature derived from single-cell RNA sequencing.

Multiple myeloma (MM) is a heterogenous plasma cell malignancy, for which the established prognostic models exhibit limitations in capturing the full spectrum of outcome variability. Leveraging single-cell RNA-sequencing data, we developed a novel plasma cell gene signature. We evaluated and validated the associations of the resulting plasma cell malignancy (PBM) score with disease state, progression and clinical outcomes using data from five independent myeloma studies consisting of 2115 samples (1978 MM, 65 monoclonal gammopathy of undetermined significance, 35 smoldering MM, and 37 healthy controls). Overall, a higher PBM score was significantly associated with a more advanced stage within the spectrum of plasma cell dyscrasias (all p < 0.05) and a shorter overall survival in MM (hazard ratio, HR = 1.72; p < 0.001). Notably, the prognostic effect of the PBM score was independent of the International Staging System (ISS) and Revised ISS (R-ISS). The downstream analysis further linked higher PBM scores with the presence of cytogenetic abnormalities, TP53 mutations, and compositional changes in the myeloma tumor immune microenvironment. Our integrated analyses suggest the PBM score may provide an opportunity for refining risk stratification and guide decisions on therapeutic approaches to MM.

Targeting NUPR1-dependent stress granules formation to induce synthetic lethality in KrasG12D-driven tumors.

We find that NUPR1, a stress-associated intrinsically disordered protein, induced droplet formation via liquid-liquid phase separation (LLPS). NUPR1-driven LLPS was crucial for the creation of NUPR1-dependent stress granules (SGs) in pancreatic cancer cells since genetic or pharmacological inhibition by ZZW-115 of NUPR1 activity impeded SGs formation. The KrasG12D mutation induced oncogenic stress, NUPR1 overexpression, and promoted SGs development. Notably, enforced NUPR1 expression induced SGs formation independently of mutated KrasG12D. Mechanistically, KrasG12D expression strengthened sensitivity to NUPR1 inactivation, inducing cell death, activating caspase 3 and releasing LDH. Remarkably, ZZW-115-mediated SG-formation inhibition hampered the development of pancreatic intraepithelial neoplasia (PanINs) in Pdx1-cre;LSL-KrasG12D (KC) mice. ZZW-115-treatment of KC mice triggered caspase 3 activation, DNA fragmentation, and formation of the apoptotic bodies, leading to cell death, specifically in KrasG12D-expressing cells. We further demonstrated that, in developed PanINs, short-term ZZW-115 treatment prevented NUPR1-associated SGs presence. Lastly, a four-week ZZW-115 treatment significantly reduced the number and size of PanINs in KC mice. This study proposes that targeting NUPR1-dependent SGs formation could be a therapeutic approach to induce cell death in KrasG12D-dependent tumors.

Verteporfin induces lipid peroxidation and ferroptosis in pancreatic cancer cells.

Pancreatic ductal adenocarcinoma (PDAC) has extremely poor prognosis, with a 5-year survival rate of approximately 11 %. Yes-associated protein (YAP) is a major downstream effector of the Hippo-YAP pathway and plays a pivotal role in regulation of cell proliferation and organ regeneration and tumorigenesis. Activation of YAP signaling has been associated with PDAC progression and drug resistance. Verteporfin (VP) is a photosensitizer used for photodynamic therapy and previous work showed that it can function as a YAP inhibitor. The efficacy of VP on human cancer are being tested in several trials. In this study, we examined the effect of VP on reactive oxygen species (ROS) and lipid peroxidation in pancreatic cancer cells, by using fluorescent molecular probes and by measuring the levels of malondialdehyde, a metabolic byproduct and marker of lipid peroxidation. We found that VP causes rapid increase of both overall ROS and lipid peroxide levels, independent of light activation. These effects were not dependent on YAP, as knockdown of YAP did not cause ROS or lipid peroxidation or enhance VP-induced ROS production. Temoporfin, another photodynamic drug, did not show similar activities. In addition, VP treatment led to loss of cell membrane integrity and reduction of viability. Notably, the activity of VP to induce lipid peroxidation was neutralized by ferroptosis inhibitors ferrostatin-1 or liproxstatin-1. VP treatment also reduced the levels of glutathione peroxidase 4 (GPX4), an enzyme that protects against lipid peroxidation. These results indicate that VP can induce lipid peroxidation and ferroptosis in the absence of light activation. Our findings reveal a novel mechanism by which VP inhibits tumor growth and provide insights into development of new therapeutic strategies for the treatment of pancreatic cancer.

Deep computational phenotyping of genomic variants impacting the SET domain of KMT2C reveal molecular mechanisms for their dysfunction.

Introduction: Kleefstra Syndrome type 2 (KLEFS-2) is a genetic, neurodevelopmental disorder characterized by intellectual disability, infantile hypotonia, severe expressive language delay, and characteristic facial appearance, with a spectrum of other distinct clinical manifestations. Pathogenic mutations in the epigenetic modifier type 2 lysine methyltransferase KMT2C have been identified to be causative in KLEFS-2 individuals. Methods: This work reports a translational genomic study that applies a multidimensional computational approach for deep variant phenotyping, combining conventional genomic analyses, advanced protein bioinformatics, computational biophysics, biochemistry, and biostatistics-based modeling. We use standard variant annotation, paralog annotation analyses, molecular mechanics, and molecular dynamics simulations to evaluate damaging scores and provide potential mechanisms underlying KMT2C variant dysfunction. Results: We integrated data derived from the structure and dynamics of KMT2C to classify variants into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). When compared with controls, these variants show values reflecting alterations in molecular fitness in both structure and dynamics. Discussion: We demonstrate that our 3D models for KMT2C variants suggest distinct mechanisms that lead to their imbalance and are not predictable from sequence alone. Thus, the missense variants studied here cause destabilizing effects on KMT2C function by different biophysical and biochemical mechanisms which we adeptly describe. This new knowledge extends our understanding of how variations in the KMT2C gene cause the dysfunction of its methyltransferase enzyme product, thereby bearing significant biomedical relevance for carriers of KLEFS2-associated genomic mutations.

An EMR-Based Approach to Determine Frequency, Prescribing Pattern, and Characteristics of Patients Receiving Drugs with Pharmacogenomic Guidelines.

(1) Background: This retrospective analysis utilizing electronic medical record (EMR) data from a tertiary integrated health system sought to identify patients and prescribers who would benefit from pharmacogenomic (PGx) testing based on Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines. (2) Methods: EMR data from a clinical research data warehouse were analyzed from 845,518 patients that had an encounter between 2015 and 2019 at an academic medical center. Data were collected for 42 commercially available drugs with 52 evidence-based PGx guidelines from CPIC. Provider data were obtained through the EMR linked by specialty via national provider identification (NPI) number. (3) Results: A total of 845,518 patients had an encounter in the extraction period with 590,526 medication orders processed. A total of 335,849 (56.9%) patients had medication orders represented by CPIC drugs prescribed by 2803 providers, representing 239 distinct medications. (4) Conclusions: The results from this study show that over half of patients were prescribed a CPIC actionable medication from a variety of prescriber specialties. Understanding the magnitude of patients that may benefit from PGx testing, will enable the development of preemptive testing processes, physician support strategies, and pharmacist workflows to optimize outcomes should a PGx service be implemented.

A multi-layered computational structural genomics approach enhances domain-specific interpretation of Kleefstra syndrome variants in EHMT1.

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretations of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within the domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism for variant dysfunction.

RAG genomic variation causes autoimmune diseases through specific structure-based mechanisms of enzyme dysregulation.

Interpreting genetic changes observed in individual patients is a critical challenge. The array of immune deficiency syndromes is typically caused by genetic variation unique to individuals. Therefore, new approaches are needed to interpret functional variation and accelerate genomics interpretation. We constructed the first full-length structural model of human RAG recombinase across four functional states of the recombination process. We functionally tested 182 clinically observed RAG missense mutations. These experiments revealed dysfunction due to recombinase dysfunction and altered chromatin interactions. Structural modeling identified mechanical and energetic roles for each mutation. We built regression models for RAG1 (R2 = 0.91) and RAG2 (R2 = 0.97) to predict RAG activity changes. We applied our model to 711 additional RAG variants observed in population studies and identified a subset that may impair RAG function. Thus, we demonstrated a fundamental advance in the mechanistic interpretation of human genetic variations spanning from rare and undiagnosed diseases to population health.

Writers and readers of H3K9me2 form distinct protein networks during the cell cycle that include candidates for H3K9 mimicry.

Histone H3 lysine 9 methylation (H3K9me), which is written by the Euchromatic Histone Lysine Methyltransferases EHMT1 and EHMT2 and read by the heterochromatin protein 1 (HP1) chromobox (CBX) protein family, is dysregulated in many types of cancers. Approaches to inhibit regulators of this pathway are currently being evaluated for therapeutic purposes. Thus, knowledge of the complexes supporting the function of these writers and readers during the process of cell proliferation is critical for our understanding of their role in carcinogenesis. Here, we immunopurified each of these proteins and used mass spectrometry to define their associated non-histone proteins, individually and at two different phases of the cell cycle, namely G1/S and G2/M. Our findings identify novel binding proteins for these writers and readers, as well as corroborate known interactors, to show the formation of distinct protein complex networks in a cell cycle phase-specific manner. Furthermore, there is an organizational switch between cell cycle phases for interactions among specific writer-reader pairs. Through a multi-tiered bioinformatics-based approach, we reveal that many interacting proteins exhibit histone mimicry, based on an H3K9-like linear motif. Gene ontology analyses, pathway enrichment, and network reconstruction inferred that these comprehensive EHMT and CBX-associated interacting protein networks participate in various functions, including transcription, DNA repair, splicing, and membrane disassembly. Combined, our data reveals novel complexes that provide insight into key functions of cell cycle-associated epigenomic processes that are highly relevant for better understanding these chromatin-modifying proteins during cell cycle and carcinogenesis.

A Multi-Layered Computational Structural Genomics Approach Enhances Domain-Specific Interpretation of Kleefstra Syndrome Variants in EHMT1.

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretation of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within this domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism of variant dysfunction.

Beyond structural bioinformatics for genomics with dynamics characterization of an expanded KRAS mutational landscape.

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations, which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

Advanced computational analysis of CD40LG variants in atypical X-linked hyper-IgM syndrome.

X-Linked Hyper-IgM Syndrome is caused by pathogenic variants in CD40LG. Three patients with atypical clinical and immunological features were identified with variants in CD40LG requiring further characterization. Flow cytometry was used to evaluate CD40L protein expression and binding capacity to a surrogate receptor, CD40-muIg. Though functional anomalies were observed, there was still a lack of clarity regarding the underlying mechanism. We developed structural models for wild-type and the three variants of CD40L protein observed in these patients (p. Lys143Asn, Leu225Ser and Met36Arg) to evaluate structural alterations by molecular mechanic calculations, and assess protein movement by molecular dynamic simulations. These studies demonstrate that functional analysis of variants of unknown significance in CD40LG can be supplemented by advanced computational analysis in atypical clinical contexts. These studies in combination identify the deleterious effects of these variants and potential mechanisms for protein dysfunction.

KLF5 and p53 comprise an incoherent feed-forward loop directing cell-fate decisions following stress.

In response to stress, cells make a critical decision to arrest or undergo apoptosis, mediated in large part by the tumor suppressor p53. Yet the mechanisms of these cell fate decisions remain largely unknown, particularly in normal cells. Here, we define an incoherent feed-forward loop in non-transformed human squamous epithelial cells involving p53 and the zinc-finger transcription factor KLF5 that dictates responses to differing levels of cellular stress from UV irradiation or oxidative stress. In normal unstressed human squamous epithelial cells, KLF5 complexes with SIN3A and HDAC2 repress TP53, allowing cells to proliferate. With moderate stress, this complex is disrupted, and TP53 is induced; KLF5 then acts as a molecular switch for p53 function by transactivating AKT1 and AKT3, which direct cells toward survival. By contrast, severe stress results in KLF5 loss, such that AKT1 and AKT3 are not induced, and cells preferentially undergo apoptosis. Thus, in human squamous epithelial cells, KLF5 gates the response to UV or oxidative stress to determine the p53 output of growth arrest or apoptosis.

Priming therapy by targeting enhancer-initiated pathways in patient-derived pancreatic cancer cells.

BACKGROUND: Systems biology leveraging multi-OMICs technologies, is rapidly advancing development of precision therapies and matching patients to targeted therapies, leading to improved responses. A new pillar of precision oncology lies in the power of chemogenomics to discover drugs that sensitizes malignant cells to other therapies. Here, we test a chemogenomic approach using epigenomic inhibitors (epidrugs) to reset patterns of gene expression driving the malignant behavior of pancreatic tumors. METHODS: We tested a targeted library of ten epidrugs targeting regulators of enhancers and super-enhancers on reprogramming gene expression networks in seventeen patient-derived primary pancreatic cancer cell cultures (PDPCCs), of both basal and classical subtypes. We subsequently evaluated the ability of these epidrugs to sensitize pancreatic cancer cells to five chemotherapeutic drugs that are clinically used for this malignancy. FINDINGS: To comprehend the impact of epidrug priming at the molecular level, we evaluated the effect of each epidrugs at the transcriptomic level of PDPCCs. The activating epidrugs showed a higher number of upregulated genes than the repressive epidrugs (χ2 test p-value <0.01). Furthermore, we developed a classifier using the baseline transcriptome of epidrug-primed-chemosensitized PDPCCs to predict the best epidrug-priming regime to a given chemotherapy. Six signatures with a significant association with the chemosensitization centroid (R ≤ -0.80; p-value < 0.01) were identified and validated in a subset of PDPCCs. INTERPRETATION: We conclude that targeting enhancer-initiated pathways in patient-derived primary cells, represents a promising approach for developing new therapies for human pancreatic cancer. FUNDING: This work was supported by INCa (Grants number 2018-078 to ND and 2018- 079 to JI), Canceropole PACA (ND), Amidex Foundation (ND), and INSERM (JI).

Beyond Structural Bioinformatics for Genomics with Dynamics Characterization of an Expanded KRAS Mutational Landscape.

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

Prognostic landscape of mitochondrial genome in myelodysplastic syndrome after stem-cell transplantation.

Despite mitochondrial DNA (mtDNA) mutations are common events in cancer, their global frequency and clinical impact have not been comprehensively characterized in patients with myelodysplastic neoplasia (also known as myelodysplastic syndromes, MDS). Here we performed whole-genome sequencing (WGS) on samples obtained before allogenic hematopoietic cell transplantation (allo-HCT) from 494 patients with MDS who were enrolled in the Center for International Blood and Marrow Transplant Research. We evaluated the impact of mtDNA mutations on transplantation outcomes, including overall survival (OS), relapse, relapse-free survival (RFS), and transplant-related mortality (TRM). A random survival forest algorithm was applied to evaluate the prognostic performance of models that include mtDNA mutations alone and combined with MDS- and HCT-related clinical factors. A total of 2666 mtDNA mutations were identified, including 411 potential pathogenic variants. We found that overall, an increased number of mtDNA mutations was associated with inferior transplantation outcomes. Mutations in several frequently mutated mtDNA genes (e.g., MT-CYB and MT-ND5) were identified as independent predictors of OS, RFS, relapse and/or TRM after allo-HCT. Integration of mtDNA mutations into the models based on the Revised International Prognostic Scores (IPSS-R) and clinical factors related to MDS and allo-HCT could capture more prognostic information and significantly improve the prognostic stratification efforts. Our study represents the first WGS effort in MDS receiving allo-HCT and shows that there may be clinical utility of mtDNA variants to predict allo-HCT outcomes in combination with more standard clinical parameters.